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INTRODUCTION 

A system A with energy E and degeneracy 

PA{E)c^tx^{kAE) 
while seemingly having a partition function of the form 

Z(r)= / pAiE)exp{-E/T)dE 



(1) 



(2) 

for all temperatures T < in fact admits only one temperature T = Ta = 1/kA and it imparts that temperature to 
any system coupled to it. 

The partition function of Eq. (|2j implies that an external thermostat/? which, by definition has Pb{E) °^ exp(— ^s/i), 
can impart its temperature Tb — I / ks to the system A. This is not so, as can be seen by considering the generating 
micro-canonical partition 



The most probable partition is given by 
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But this is hardly possible since in general Ta^Tb'- two thermostats can never be at equilibrium unless they are at the 
same temperature. 

This preamble is motivated by the fact that the empirical hadronic mass spectra (Hagedorn spectra fP, '5]), the 
Statistical Bootstrap Model (SBM) |3, 4, 5] and the MIT bag model |6] have a degeneracy whose leading term is 
of the form of Eq. Q. It is the aim of this paper to explore in a pedagogical manner the implications of such a 
spectrum, making only passing references to the more complex physical situations occurring in particle-particle and 
nucleus-nucleus collisions. 

Hagedorn noted that the hadronic mass spectrum (level density) has the asymptotic (m ^ o°) form 



P^(ot) wexp(m/r/^) 



(5) 



where m is the mass of the hadron in question and is the temperature associated with the mass spectrum (Hl^]- 
The question of the mass range over which (|5} is vaHd is still under discussion 

The M.I.T bag model f6'| of partonic matter reproduces this behavior via a constant pressure B of a "bag" of partonic 
matter 00]. The pressure p inside a bag at equilibrium without additional conserved quantities is 



P=W'^s-B, (6) 

where g is the number of partonic degrees of freedom. The bag constant forces a constant temperature Tb on the bag. 
Similarly, the enthalpy density e of the bag 

y 30 ^ 

is constant. Here H is the enthalpy and V is the volume of the bag. Thus, an injection of an arbitrary amount of energy 
leads to an isothermal, isobaric expansion of the bag and the bag entropy S is proportional to H: 

T Jo T Tb' 

where 52 is the change in heat of the bag. The bag's spectrum (level density) is then p — exp (5') given by Eq. (jSj with 

Tb = T.ye and H = m. 

Following our recent results f9',' 10], we show here that a system possessing a Hagedorn-like spectrum, 
characterized by an entropy of the form (|8j, not only has a unique microcanonical temperature 
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Tb, (9) 



but also imparts this same temperature to any other system to which Jf is coupled. In the language of standard 
thermodynamics: is a perfect thermostat. 

The property of a perfect thermostat is well known. For instance, it is indifferent to the transfer of any portion of 
its energy to any parcel within itself, no matter how small. In other words, it is at the limit of phase stability and its 
internal fluctuations of the energy density are maximal. 



HARMONIC OSCILLATOR COUPLED TO 

In order to demonstrate the thermostatic behavior of a Hagedorn system, let us begin by coupling to a one 
dimensional harmonic oscillator and use a microcanonical treatment. The unnormalized probability P{e) for finding 
an excitation energy e in the harmonic oscillator out of the system's total energy E is 

P(e) - p._^(£-e)posc(e) 

= exp(^^^ =p,^(£)exp(^-^^. (10) 

Recall that for a one dimensional harmonic oscillator Posc is a constant. The energy spectrum of the oscillator is 
canonical up to the upper limit Emax — E with an inverse slope (temperature) of independent of E. The mean value 
of the energy of the oscillator is given by 



£^T^ 



exp{E/T^)-l 



(11) 



Thus in the limit that E ^ e ^ T^, i.e. no temperature other that Tjf is admitted. In the standard language 
of statistical mechanics this example means that a one dimensional harmonic oscillator can be used as an ideal 
thermometer. 



AN IDEAL VAPOR COUPLED TO M' 



For a more physically relevant example, let us consider a vapor of ^ 1 non-interacting Boltzmann particles of mass 
niB and degeneracy gB coupled to M' . The microcanonical level density of the vapor with kinetic energy e is 

where V is is the volume. The microcanonical partition of the total system is 



Ptotal (£■,£) = P^(£'-e)Pvapor(e) 



(13) 



Just as with the harmonic oscillator, the distribution of the vapor is exactly canonical up to — E, if the particles 
are independently present, or e„^ax = E — mN, if the particles are generated by Jf. In either case, the temperature of 
the vapor is always T/^. 

The maximum of ptotai {E, e) with respect to e gives the most probable kinetic energy per particle as 

dp,ME,e)^3N_^^ e_^l 

de 2e Tyr N 2^' ^ ' 

provided that E > msN + jNT^^. (For mbN < E < mbN + jNT^, the most probable value of the kinetic energy 
per particle is ^ = ^ — < jT^; for E < mbN, = 0. ) Again is the sole temperature characterizing the 
distribution up to the microcanonical cut-off, which may be above or below the maximum of the distribution since the 
form of ptotai(£^, e) is independent of E. 

The maximum of Ptotai(£, e) with respect to at fixed V is given by 
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where Eq. MAI was used for £. Thus the most probable particle density of the vapor is independent of V: 



gB\— e ^x- =n,^. (16) 



y ° \ 1% ^ 

Equation (I16> demonstrates that not only is M' a perfect thermostat but also a perfect particle reservoir. Particles 
of different mass m will be in chemical equilibrium with each other At equilibrium, particles are emitted from 
and form a saturated vapor at coexistence with M' at temperature T^f. This describes a first order phase transition 
(hadronic to partonic). Coexistence occurs at a single temperature fixed by the bag pressure. 

These results explain the common value of: the hadronization temperatures obtained within the statistical hadroniza- 
tion model 1 1 1 1; the inverse slopes of the transverse mass spectra of hadrons observed in high energy elementary par- 
ticle collisions 1 12, 13|; and the transition temperature from lattice QCD calculations for low baryonic density 
For further discussion see [10]. 



^ AS A RADIANT BAG 

Let us assume that M' is a bag thick enough to absorb any given particle of the vapor striking it. Then, detailed balance 
requires that on average M' radiates back the same particle. Under these conditions particles can be considered to be 
effectively emitted from the surface of M' . Thus the relevant fluxes do not depend in any way upon the inner structure 
of Jf". 

In fact, the results given in equations ( I14> and (I16> show that the saturated vapor concentration depends only upon 
niB and as long as ^ is present. A decrease in the volume V does not increase the vapor concentration, but 



induces a condensation of the corresponding amount of energy out of the vapor and into . An increase in V keeps 
the vapor concentration constant via evaporation of the corresponding amount of energy out of M' and into the vapor 
This is reminiscent of liquid-vapor equihbrium at fixed temperature, except that here coexistence occurs at a single 
temperature , rather than over a range of temperatures as in ordinary fluids. 

The bag wall is Janus faced: one side faces the partonic world, and, aside from conserved charges, radiates a partonic 
black body radiation responsible for balancing the bag pressure; the other side faces the hadronic world and radiates 
a hadronic black body radiation, mostly pions. Both sides of the bag wall are at the temperature T,^. It is tempting 
to attribute most, if not all, of the hadronic and partonic properties to the wall itself, possibly even the capability to 
enforce conservation laws globally (quantum number conductivity). Despite the fact that this wall is an insurmountable 
horizon, with hadronic measurements such as bag size and total radiance we can infer some properties of the partonic 
world, e.g. the number of degrees of freedom 1 12]. 

We can estimate an upper limit for the emission time using the outward energy flux of particles radiated from the 
bag. At equihbrium the in-going and out-going fluxes must be the same, thus the outward flux of particles in the 
nonrelativistic approximation using Eq. J16t is 



Using the technique developed in one finds the energy flux ^>E.yf and momentum flux p-^^A as 

(^Eyf ^ + 2 ) (^nyy , /7rad = ^ n.ir 7> • (18) 

The pressure /Jrad exerted on the bag by its radiation can be compared to the intrinsic bag pressure in Eq. (|6j: for pions 
Prad ~ 0.02B. The time T for the bag to dissolve into its own radiation is approximately 

(19) 
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where is the initial bag radius and £0 is the initial bag total energy. 

The fluxes written in Eqs. (^} and ( I18> (particle or energy per unit surface area) are integrated over an assumed 
spherical bag to give the result in Eq. J19> . However, because of the lack of surface tension, the bag's maximum 
entropy corresponds to either an elongated (cylinder) or a flattened shape (disc). Thus, Eq. ( I19> should be interpreted 
as an upper limit. More detailed studies of hadron emission from bags concerning hydrodynamic shock waves and 
freeze out shocks can be found elsewhere 

The decoupling between the vapor concentration and and Tyf occurs when M' has completely evaporated (i.e. 



when E — niBN— jNTj^ = 0) at a volume of 



1 E 

Vrf- 3-—. (20) 



The disappearance of Jff allows the vapor concentration to decrease inversely proportionally to V as 

V ~ V 



(21) 



The temperature, however, remains fixed at Tjf because of conservation of energy and particle number above Vd- 
Solid curves in Fig. 1 show this schematically. 

The discussion above assumes that the Hagedorn spectrum extends down to m ~ 0. However, experimentally there 
appears to be a lower cut off of the spectrum at mo. This modifies the above results as follows (for a detailed analysis 
see the section "Generalization to a Complete Hagedorn Spectrum"). 

For energies E — niBN — e ^ otq and V < the above results hold as written. However, if we increase the volume 
well beyond Vd at which the Hagedorn spectrum is truncated at mo, the situation is slightly different. evaporates 
until its mass is mo- If the entire mass of J^f is fully transformed into vapor particles as the volume is increased further, 
then the excess particles temporarily increase the concentration and permanently decrease the temperature. As the 
volume increases further, the concentration changes inversely proportional to V 
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FIGURE 1. Typical behavior of the entire system's temperature T and concentration A^/V as the function of the 
system's volume V in the absence of restrictions (solid curve) and for a finite cut off at mo of the Hagedorn spectrum 
(dashed curve). 



while the temperature remains constant at 



T = 



(23) 



Dashed curves in Fig. 1 show this schematically. 



FRAGMENTATION OF M' 



A question of interest is the stability of ^ against fragmentation. If the translational degrees of freedom are neglected, 
M' is indifferent to fragmentation into an arbitrary number of particles of arbitrary mass (within the constraints of 
mass/energy conservation). 

Let us now consider the case in which the mass of the vapor particle is allowed to be free. The system's level 
density Ptotai(£^,e) is still given by Eq. (I13> . Using Eqs. ( I14> and ( I16> . one finds the most probable value of the system's 
level density as pj^(^j(£',e) w exp[5'*], where the entropy is S* = E jT,^^ +N. Differentiating Ptotai(^'^) ^^'-^ respect 
to niB and applying Eq. ( I16t gives 
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(24) 



i.e. the last equality provides the maximum of level density for N ^0. Since all the intrinsic statistical weights in 
Ptotai(^' e) are factored into a single Jf, the system breaks into fragments with mb — except for one whose mass 
is determined by mass/energy conservation. 



Substituting the most probable value of e and ms into the most probable value of one obtains the vapor 
concentration 

N ( ' 
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(25) 



The density of the vapor of nonrelativistic particles acquires the form typical of the ultrarelativistic limit. 

If the value of mass given by Eq. i24\ does not exist, then the most probable value of level density p*g^^^{E,e) 
corresponds to the mass m* which is nearest to jTjf and N{m*) given by Eq. il6l . In terms of hadron spectroscopy 
the value of m* that maximizes the level density p*g^^^{E,e) is the pion mass. 

If Jf is required to fragment totally into a number of equal fragments of mass mn all endowed with their 
translational degrees of freedom, then (for = 1) 
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(26) 



where in the last step we substituted the most probable value of the kinetic energy MA\ and used the Stirling formula 



for From Eq. (|26) it is seen that all the Hagedorn factors collapse into a single one with the m-independent 



argument E. Maximization of (I26> with respect to nin leads to 



d\npj 



3N 



dniH 2mH 



= 0, 



(27) 



which is consistent with = 1 and mn ~ E, namely a single Hagedorn particle with all the available mass. 

This again illustrates the indifference of Jff toward fragmentation. Of course Eq. (I14> gives directly the mass 
distribution of the Hagedorn fragments under the two conditions discussed above. These results justify the assumption 
of the canonical formulation of the statistical hadronization model that smaller clusters appear from a single large 
clustering. 



INTERMEDIATE CONCLUSIONS 

A system jyf, with a Hagedorn-like mass spectrum, is a perfect thermostat and a perfect particle reservoir. Conse- 
quently, any system coupled to jyf can have only the temperature of : Tjf. This behavior may explain the common 
value of: the hadronization temperatures obtained within statistical models; the transition temperature from lattice 
QCD calculations for low baryonic density; and the inverse slopes of the transverse mass spectra of hadrons (temper- 
ature) observed in high energy elementary particle collisions and high energy nucleus-nucleus collisions (for details 
see 1 10]). The common temperature of the experimental spectra suggest that the observed particles originate from an 
J^-liks system. 

The hadronic side of Jf radiates particles in preexisting physical and chemical equilibrium just as a black body 
radiates photons in physical and chemical equilibirum (compare to Ref. |20]). Particles emitted from form a 
saturated vapor that coexists with J/f. This coexistence describes a first order phase transition (hadronic to partonic) 
and occurs at a single temperature fixed by the bag pressure. An system is nearly indifferent to fragmentation into 
smaller J^f systems. A lower cut-off in the mass spectrum does not alter our results 1 101 . 



GENERALIZATION TO A COMPLETE HAGEDORN SPECTRUM 

To have a more realistic model we should consider a more complicated Hagedorn mass spectrum g//(m//) = 
exp[m///rj)f ] (mo/m//)" for the resonance masses nin above the lower cut-off nio 3> Tjf (a is a parameter discussed 
below). Let us study the microcanonical ensemble of Nb Boltzmann point-like particles of mass mg and degeneracy 
gB, and Nh hadronic point-like resonances of mass nin with a mass spectrum g//(m//) assuming that nio > niB- A 
recent analysis |21 ] suggests that the Hagedorn mass spectrum can be established for nio <2 GeV. 

In the Statistical Bootstrap Model (SBM) |22] and the MIT bag model |0] it was found that for m// oo the 
parameter a < 3. For finite resonance masses the value of a is unknown, so it will be considered as a fixed parameter. 



The microcanonical partition of the system, with volume V, total energy E and zero total momentum, can be written 
as follows 
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(28) 



where the quantity = £{mH,Qi) ye^ ~ e{mB,Pj) and e{M,P) = \/M^ +P^j denotes the energy of the Hagedorn 

(Boltzmann) particle with the 3-momentum Q, iPj)- In order to simplify the presentation of our idea, Eq. j28> accounts 
for energy conservation only and neglects momentum conservation. 

The microcanonical partition ( I28> can be evaluated by the Laplace transform in total energy E [23]. Then the 
momentum integrals in ( I28> are factorized and can be performed analytically. The inverse Laplace transform in the 
conjugate variable A can be done analytically for the nonrelativistic and ultrarelativistic approximations of the one- 
particle momentum distribution function 
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where the auxiliary integral is denoted as 
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(30) 



Since the formal steps of further evaluation are similar for both cases, we discuss in detail the nonrelativistic limit 
only, and later present the results for the other case. The nonrelativistic approximation (MRe{X) ^ 1) for Eq. ( I28> is 
as follows 
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(31) 



where Ej^i^ = E — mnNn — niBNB is the kinetic energy of the system. 

As shown below, the most realistic case corresponds to the nonrelativistic treatment of the Hagedorn resonances 
because the resulting temperature is much smaller than their masses. Therefore, it is sufficient to consider the ultra- 
relativistic limit for the Boltzmann particles only. In this case {MReiX) ^1) the equation ( I28> can be approximated 
as 
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(32) 



where the kinetic energy does not include the rest energy of the Boltzmann particles, i.e. Eiii„ = E — iuhNh- 

Within our assumptions the above results are general and can be used for any number of particles, provided 
Nh +Nb >2. It is instructive to consider first the simplest case Nh ~ 1. This formulation of the model, in which 
a Hagedorn thermostat is always present, allows us to study the problem rigorously and provides us with a qualitative 
picture for Nh > 1 . For Nh = 1 and Nb ^ 1 we treat the mass of Hagedorn thermostat mn as a free parameter and 
determine the value which maximizes the entropy of the system. The solution m]j > of the extremum condition 
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provides the maximum of the system's entropy, if for mn = the second derivative is negative 

5^\nQ.„r{NH = 1) 



HNb+I) 
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< 0. 



(33) 



(34) 



The inequaUty i34l is a necessary condition of the maximum of the microcanonical partition. Postponing the analysis 
of (|33 till the next section, where we study it in more details, let us assume for a moment that the inequality (|33 is 



satisfied. Then the extremum condition ( I33> defines the temperature of the system of {Nb + 1 ) nonrelativistic particles 



T^*/ *\ ^ Eiii„ Tj^ /OCX 

^ <'»«'= si^^TTTF^- 

Thus, as mil ^ °° it follows that T* {m*fj) T.ye, while for finite m*n ^ Tjf and a > 5 (a < |) the temperature of the 
system is slightly larger (smaller) than the Hagedorn temperature, i.e. T* > (T* < T^). Formally, the temperature 
of the system in equation \35l may differ essentially from Tj^ for a light thermostat, i.e. for m*^ < T^. However, it is 
assumed that the Hagedorn mass spectrum exists above the cut-off mass mo ^ Tjf, thus m* ^ T^. 

THE ROLE OF THE MASS CUT-OFF 

Now we study the effect of the mass cut-off of the Hagedorn spectrum on the inequality J34I I in more detail. For a < | 
the condition ( I34> is satisfied. For a > 5 the inequality (I34> is equivalent to the following inequality 

^^^j^>^«. + i)r«), ,36, 

which means that a Hagedorn thermostat should be massive compared to the kinetic energy of the system. 
A more careful analysis shows that for a negative value of the determinant D„r {N = Nb — fa) 

D„r = {E-msNB-lT^Nf- 

4(fl-|) (E-mBNB) < 0, (37) 

equation i33\ has two complex solutions, while for D,,, — there exists a single real solution of ( I33> . Solving (I37> for 
{E — mBNs), shows that for^^V > |fl — 1, i.e. for Nb > ^,a—l the inequality ( I37t does not hold andD„r > 0. Therefore, 
in what follows we will assume that Nb > — I and only analyze the case D„r > 0. For this case equation i33l has 
two real solutions 

m^^UE-mBNB-^TrfN±VD::^- (38) 



2 

For a < J only solution is positive and corresponds to a maximum of the microcanonical partition £2„, . 
For a> J both solutions of i33\ are positive, but only is a maximum. From the two limiting cases: 

- " a)-^ for mn^O, (39) 
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for £fa„«0, (40) 



5 mn ^ 
and the fact that mjj obey the inequalities 

< mff < m^ < E-mBNB, (41) 

it is clear that m^/ — is a local minimum of the microcanonical partition n„,-, while m*fj — m^ is a local maximum 
of the partition £l„r- 

Using Eq. \3%\ for m^, it is clear that for any value of a the constraint m^ > mo is equivalent to the inequality 

E - [Jp^ - a]T*{mo) 

Nb < N'b'" ^ , . (42) 

mB + 5 T*(mo) 

Thus, at fixed energy E for all Nb < N^'" at = m+ there is a local maximum of the microcanonical partition il^^- 
with the temperature T — T*{m'^). For Nb > Nb" the maximum of the partition Q.„r cannot be reached due to the 
cut-off constraint and, consequently, the most probable state corresponds to mn — mo with T <T* (mo) from Eq. ( I35t . 
In other words, for Nb > Ng" the amount of energy E is insufficient for the mass of the Hagedorn thermostat to be 
above the cut-off mo and simultaneously maintain the temperature of the Boltzmann particles according to Eq. ( I35t . 
By assumption there is a single Hagedorn thermostat in the system, therefore, as Nb grows the temperature of the 



system decreases from T* {nio) value. Thus, the equahty ( I42t defines the kinematical Hmit for reaching the maximum 
of the microcanonical partition. 

To prove that the maximum of the microcanonical partition at mh — mjj is global it is sufficient to show that 
the constraint > mo is not consistent with the condition mjj > mo- For « < | the maximum is global because 
for < niH < mjj (mn > ) the partition Q.„r{NH = l,m//) monotonically increases (decreases) with mn- For 
a > I it is clear that the maximum at mn = is local, if the state with mass mn = mo is more probable, i.e. 
^nr{NH = l,'«o) > ^nr{NH = 1, '«//)• Due to (I41> this Can occur, if > Substituting Eq. (I38> into the last 
inequality, shows that this inequality reduces to the condition Ng > Ng". This contradicts the constraint > Mq in 
the form of Eq. (I42> . Thus, the maximum of the microcanonical partition is global. 

To complete our consideration of the nonrelativistic case let us express the partition Q in terms of the temperature 
(l35l . Applying the Stirling approximation to the factorial ( | {Nb + 1 ) — 1 ) ! for N^'" >Nb^ i and reversing the integral 
representations (|55Jl and for A = 1 /r*(m^), one finds 
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This is just the grand canonical partition of (Nb + 1) Boltzmann particles with temperature T*{mjj). UNb > Ng" ^ 1, 
then T* (m+) in should be replaced by ToiNB) = ^'^3('a^%"'°' ■ 

Fig. 1 shows that for a > j the system's temperature T — T*{m^) as a function of Nb remains almost constant for 
Nb < Ng" , reaches a maximum at Ng" and rapidly decreases like T = To (Nb ) for Nb > Ng" . For a < j the temperature 
has a plateau T = T*{m^) f or Nb < Ng", and rapidly decreases for A^g > A^^'" according to To{Nb)- 

The same results are valid for the ultrarelativistic treatment of Boltzmann particles. Comparing the nonrelativistic 
and ultrarelativistics expressions for the microcanonical partition, i.e. equations Q and i32\ . respectively, one finds that 
the derivation of the ultrarelativistic limit requires only the substitution Nb ^ 2Nb and itib /T/f — > in equations ( 1331 - 
I43> . Note that this substitution does not alter the expression for the temperature of the system, i.e. the right hand side 
of (Ol}. 

Finally, we show that for a heavy Hagedorn thermostat (m^ 3> nio) these results remain valid for a single Hagedorn 
thermostat split into Nh pieces of the same mass. Substituting nifj nifjNn in the nonrelativistic expressions Q and 
minimizing it with respect to mn, the temperature of the system in the form of equation i35\ is T*{m*fjNH), where 
the mass of Nh Hagedorn thermostats is related to the solution of equation ( I38> as m*fj = m^/NH- Since the 
original single thermostat of mass was assumed to be heavy, it follows T*{m*uNH) = T*{m^) T_^'. A more 
careful study (see also |9]) using an exact expression for the microcanonical partition of Nh Hagedorn thermostats 
of the same mass mn gives the same result, if mn ^ '«o. A generalization of these statements to the case of Nh 
heavy Hagedorn thermostats of different masses also leads to the same result. Thus, splitting a single heavy Hagedorn 
thermostat into an arbitrary number of heavy resonances (heavier than mo) does not change the temperature of the 
system. 

THE BAG SURFACE 

The bag expressions reported above contain only volume terms. Given the finite size of the bags that are typically 
considered (resonances), it may be of interest to consider finite size effects and their role in the description of the bags 
properties. The simplest generalization, assuming that the bags are leptodermous (which is supported by the short 
range of hadron-hadron interaction and by the saturating properties implicit in Eqs. (|6j and Q), is the introduction of 

surface energy. This can be done phenomenologically by introducing a V 5 term in the free energy. Then the pressure 
of a spherical bag can be written as 

O A 2 , , 1 (T A 2 aJT) 

P^-T^-B--a.iT)V-. = -T^-B--^, (44) 

where fl.s(r) is the temperature dependent surface energy coefficient, R is the bag radius and a = Using the 

thermodynamic identities for the free energy F and entropy S 
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FIGURE 2. A typical behavior of the system's temperature as the function of the number of Boltzmann particles 
Nb for a = 3 and a = for the same value of the total energy E = 30 mg. Due to the thermostatic properties of a 
Hagedorn resonance the system's temperature is nearly constant up to the kinematically allowed value A^|"' given 
by (|43. 



one can find all thermodynamic functions as follows 
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(46) 
(47) 
(48) 



In evaluating the expression (Hi we fixed the integration constant (an arbitrary function of T) to zero because the free 
energy should vanish for the bag of zero volume. 

While the magnitude of fls(r) is unknown, there are surprising consequences for fls(r) > 0. In Eq. J44> the surface 
term appears as an additional pressure to the bag pressure. Therefore, for a bag in a vacuum the total pressure should 
be zero, i.e. p = 0, and, consequently, the bag temperature acquires volume dependence: 



T(R) 



a 



2a,{T) 
3aR 



(49) 



When R is large we recover the previous bag temperature and the associated physics. When R becomes small, however, 
the bag temperature increases! The implications of this dependence are strange indeed. The first is the peculiar behavior 
of the bag's heat capacity. The second is the stability of the gas of bags (or lack thereof). The third is the signature of 
a bag's decay. 




FIGURE 3. A schematic volume dependence of the bag temperature (left panel), bag energy (middle panel) and its 
heat capacity (right panel) for the temperature independent surface tension (T) = > 0. The left and right panels 
show the volume dependence of the right hand side of Eqs. ( I49l l and ( I48> . respectively. The resulting heat capacity 
of the bag is negative (right panel). The parameter To is defined by the bag constant as follows: To = [3B/aY^^ . 



HEAT CAPACITY 



In the standard bag model the heat capacity is infinite: no matter how much energy is fed to the bag, its temperature 
remains constant llgl lioll . The only effect is to make the bag larger. This is completely consistent with what we observe 
in isobaric phase transitions in ordinary matter. Here the isobaric condition is produced by the bag constant, and the 
phase transition is from hadronic to partonic phase. 

Including surface effects, shows that the more energy is put into the bag, the lower its temperature becomes: i.e. 
the bag's heat capacity is negative. To illustrate how the negative heat capacity of the bag appears, let us consider a 
temperature independent surface tension: a.s(r) ~ Qq > 0. For this case, Eq. ( I49> shows that the bag temperature is 
decreasing function of its volume, whereas, according to Eq. (I48> . the energy of the bag is an increasing function of the 
bag volume. Therefore, the bag's heat capacity, defined as dE/dT, is negative. This is shown schematically in Fig. 3. 

For a formal analysis of the heat capacity of the bag it is necessary to use Eqs. j45> and (^). From these equations 
one can find the heat capacity of the bag at constant pressure Cp and at constant volume Cy as: 
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(50) 
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In evaluating the expression for Cp we used an explicit form of the derivative 

-1 



dV_ 
dT 



dp_ 
dT 



3V^ 
2a,{T) 



dp 



V 



O 2 d 
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(52) 



From Eqs. (15 0> and J51> it is clearly seen that for any T where flj (r) > there may exist a range of parameters 
for which the heat capacity Cp, corresponding to the bag equilibrium in vacuum, is negative. This leads to a "convex 
intruder " in the entropy or an unusual behavior of its second derivative: 



1 



(53) 




FIGURE 4. Left panel: The bag temperature Th as a graphical solution of Eq. (I55> for the linear T dependence of 
the bag surface tension. The left hand side of Eq. ( I55> is shown by a bi-quadratic parabola aT^ and its right hand 
side is depicted by the straight lines for different values of the bag volume V . The solution of Eq. \55\ is found as 
an intersection point between the parabola and the straight line. 

Right panel: Shows schematically the range of available temperatures of the bag for T independent (red curve) 
and for the linear T dependent (blue curve) surface tension of the bag. See text for the details. 



which becomes positive for this range of parameters. 

In the literature on this subject it is argued |39, 40, 41] that all small systems (comparable in size with the range of 
the prevailing force) should show this effect. However, we stress that a convex intruder in the bag model with surface 
tension exists not for small systems, but for large ones and does not disappear in thermodynamic limit. This behavior 
can be verified by examining the decay products of heavy resonances: heavier resonances should decay into light 
hadrons of a lower temperature (but never lower than T{R = °o)). 

Let us now demonstrate the appearance of a convex intruder in a few simple cases. First we consider the case of 

constant surface tension, i.e. fls(r) = > 0, in more detail. Substituting Aq into Eq. (I52> . one obtains that (^^^ < 0. 
Since the heat capacity at /:> = is defined as Cp = (^ff- j ^ (^) )' opposite to the sign of the derivative 



(If) ' which can be found from the expression for the energy of the bag: 



'p=0 

2 ( dE\ 2 an 

^4BV + 3aoV^ ^ — =4B + ^>0. (54) 

p=0 V 0*^/^=0 V7 

Thus, in the case of a constant surface tension the heat capacity at = is negative which corresponds to a convex 
intruder. 

Now we consider a surface tension with a linear T dependence in a spirit of the Fisher droplet mode l I33| or 
using a more elaborate approach of the recently solved "Hills and Dales Model" for surface deformations I34ll3jl : 
OsiT) = Co ^"^''y "^^ , which is defined for the temperatures not above the critical temperature T^. 

Introducing the notation B = jT^, one can rewrite the equilibrium condition of the bag p{Th) = as follows: 

aT^ = aT^ + 2cJ^—^^, (55) 

which should be solved for the bag temperature 7//(y). 

For positive values of Th the left hand side of Eq. (I55> is a monotonically increasing function of Th, whereas its 
right hand side is a monotonically decreasing function of Th (see the left panel of Fig. 4). Therefore, there can exist a 
single intersection point of these two functions for any positive value of bag volume V . Using Eq. \55\ . one can show 
that the inequality Th < 71- is always fulfilled, if Tc > Tq. Moreover, one can also show that the allowed interval of 



the bag temperatures is between To and with limiting cases Th{V ^ 0) ^ Tc and Th{V °°) ^ To (see the right 
panel of Fig. 4). Similarly from Eq. ( I52t one finds that the bag temperature decreases, while bag volume grows, i.e. 

(It) ^ ^'^^ '^^y ^ '^^^ '^^y — ' 'Siii'^s the range of allowed bag temperatures is bound between To and T^, then 
from Eq. i50\ one can immediately see that for any Th < Tc the heat capacity of the bag at /9 = is negative for large 
volumes. Thus, in the case of a linear T dependence of the surface tension of the bag the convex intruder exists for 
large volumes of the bag. In fact, this proves the following statement: if the surface tension as{T) > is a regular 

is finite provided that the solution TiV) ofEq. H491 does not vanish in the limit 



function ofT that ipfr < and 



y oo, then in this limit the heat capacity at constant p — is negative and sign Cp = sign < 0. 

STABILITY OF A GAS OF BAGS 

A gas of resonances (bags) is frequently considered either in equilibrium or in transport problems. In our previous 
papers sua (see also preceding sections) we have shown that an ordinary bag (no surface energy) is nearly indifferent 
to fragmentation into smaller bags. In fact, under rather general conditions it appears that there is a mild tendency for 
a gas of bags to collapse into a single one. We show now that the introduction of the surface leads to an even stranger 
tendency for a gas of bags toward collapse. 

Let us assume an arbitrary mass distribution in a gas of bags, and for simplicity, let us assume that the gas is confined 
in a fixed volume with its decay products (say pions). The gas cannot be isothermal since the smaller bags have larger 
temperature than the big ones. Thus the smaller bags evaporate first and their evaporation products are absorbed by the 
larger bags until only one remains. It may be argued that isothermicity can be achieved by having all the bags to be of 
the same size. But this situation is clearly unstable. Any small perturbation in size will lead to a catastrophic collapse 
of all bags into a single one. 



DECAY OF A BAG 

A hot bag, unless constrained by conserved quantities, must decay. As it decays, the instantaneous spectrum of the 
decay products indicates the bag's instantaneous temperature. Without surface effects the bag temperature is constant 
and the overall spectrum and the instantaneous spectrum is the same. 

With the surface effects, as the bag decays and becomes smaller, its temperature increases. Therefore the overall 
spectrum integrated over the overall decay must differ from the instantaneous spectrum associated with each tempera- 
ture. The shape deviation of the overall spectrum from that of an instantaneous spectrum at fixed temperature may be 
an interesting observable to characterize both the effect and the magnitude of the surface energy and its temperature 
dependence. It is amusing to notice the similarities with a black hole and its temperature as it decays through the 
Hawking radiation. 



CONCLUSIONS 

In Refs. fg", we generalized the SBM results |22] to systems of finite energy by showing explicitly that even a 
single resonance with the Hagedorn mass spectrum degeneracy, i.e. a Hagedorn thermostat, keeps an almost constant 
temperature close to Tj^f for any number of Boltzmann particles 3 < Nb < Ng" . For the high energy limit E ^ nio 
this means that a single Hagedorn resonance defines the temperature of the system to be only slightly different from 
Tyf until the energy of the Hagedorn thermostat is almost negligible compared to E. In contrast to the grand canonical 
formulation of the original SBM liSSl . in the presence of a Hagedorn thermostat the temperature Tjjf can be reached at 
any energy density. 

The thermostatic nature of a Hagedorn system obviously explains the ubiquity of both the inverse slopes of measured 
transverse mass spectra |13] and hadronization temperature found in numerical simulations of hadrons created in 
elementary particle collisions at high energies 11^ 111 in. By a direct evaluation of the microcanonical partition 
we showed that in the presence of a single Hagedorn thermostat the energy spectra of particles become exponential 
with no additional assumptions, e.g. phase space dominance |25i] or string tension fiuctuations |26|. Also the limiting 
temperature found in the URQMD calculations made in a finite box (^T^l can be explained by the effect of the Hagedorn 
thermostat. We expect that, if the string parametrization of the URQMD in a box iEtIi was done microcanonically 
instead of grand canonically, the same behavior would be found. 

The Hagedorn thermostat model generalizes the statistical hadronization model which successfully describes the 
particle multiplicities in nucleus-nucleus and elementary collisions The statistical hadronization model 



accounts for the decay of heavy resonances (clusters in terms of Refs. fl9lfTlll24ll ') only and does not consider the 
additional particles, e.g. light hadrons, free quarks and gluons, or other heavy resonances. As we showed, the splitting 
of a single heavy Hagedorn resonance into several does not change the temperature of the system. This finding justifies 
the main assumption of the canonical formulation of the statistical hadronization model 1 19] that smaller clusters may 
be reduced to a single large cluster Also our approach naturally explains why a sophisticated transport model |28], 
which treats the hadronic reactions microscopically, leads to the thermal equilibration at the Hagedorn temperature Tj^ 
and to a chemical composition of hadrons given by the equilibrium values of particle concentrations. Thus, recalling 
the MIT Bag model interpretation of the Hagedorn mass spectrum ||7|,|8|], we conclude that quark gluon matter confined 
in heavy resonances (hadronic bags) controls the temperature of surrounding particles close to T^-, and, therefore, this 
temperature can be considered as a coexistence temperature for confined color matter and hadrons. Moreover, as we 
showed, the emergence of a coexistence temperature does not require the actual deconfinement of the color degrees 
of freedom, which, in terms of the Gas of Bags Model L29.I . is equivalent to the formation of the infinitely large and 
infinitely heavy hadronic bag. 

Within the framework of the Hagedorn thermostat model we found that even for a single Hagedorn thermostat and 
a > \ the system's temperature T = r*(m^) as a function of A^b remains almost constant for A^b < A^g", reaches a 
maximum atNg" and rapidly decreases for A^g >Ng^ (see Fig. 1). For a < | the temperature has a plateau T = T*{m'^) 
for A^B < Ng", and rapidly decreases for A^b > N^" . If such characteristic behavior of the hadronization temperature or 
the hadronic inverse slopes can be measured as a function of event multiplicity, it may be possible to experimentally 
determine the value of a. For quantitative predictions it is necessary to include more hadronic species into the model, 
but this will not change our result. 

If we apply the Hagedorn thermostat model to elementary particle collisions at high energy, then, as shown above, 
the temperature of created particles will be defined by the most probable mass of the Hagedorn thermostat. If the 
most probable resonance mass grows with the energy of collision, then the hadronization temperature should decrease 
(increase) to T,^ for a > | (a < |). Such a decrease is observed in reactions of elementary particles at high energies, 
see Table 1 in Ref. Q. 

Further we discussed the effects of the surface energy on the properties of a bag (heavy resonance) in vacuum. We 
showed that in the presence of non-zero surface tension the temperature of the bag (and any system in thermal contact 
it) acquires a volume dependence, so that smaller bags are hotter. The temperature of large bags cannot be smaller than 
the Hagedron temperature. Under not too restrictive conditions we found that the heat capacity of large bags at zero 
pressure is negative, i.e. such bags have abnormal behavior of the second derivative of entropy with respect to energy. 
These unusual properties lead to an instability of any number of bags other than one. We argued that the temperature 
of the decay products of the evaporating bag should grow during the evaporation process, which, hopefully, can be 
observed. 

In order to apply these results in a more physical fashion to the quark gluon plasma formation in relativistic nucleus- 
nucleus collisions (where the excluded volume effects are known to be important |30, 29, 31, 32, 36, 37] for all 
hadrons) the eigen volumes of all particles should be incorporated into the model. For pions this should be done in 
relativistic fashion llssil . Also the effect of finite width of Hagedorn resonances may be important fs*] and should be 
studied. 
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